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Abstract: The problem of the vacuum energy decay is studied for both signs of 

^ the cosmological constant, through the analysis of the vacuum survival amphtude, 

>■ defined in terms of the conformal time, z, by A{z,z') = (vac 2 1 vac 2'). Transition 

^^ amphtudes are computed for finite time-span, Z = z' — z, and their late time behavior 

^~~* (directly related to the putative decay width of the state) as well as the transients 

• are discussed up to first order in the coupling constant, A. 
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1. Introduction 

It has been claimed [1] that the free energy corresponding to an interacting theory 
in de Sitter space has got an imaginary part that can be interpreted as some sort of 
instabihty. 

The way this imaginary part has to be (perhaps naively) computed is by doing 
the path integral at imaginary time. In Poincare coordinates of de Sitter space (dSn) 
in which the metric reads 



ds' 



z^ 



this leads to the metric of euclidean (anti) de Sitter space (EAdSn) with metric 



ds' 



^2 



(and not to the metric on the sphere 5'„), and the corresponding free energy has been 
computed by us up to the one loop order in the first paper of [2] , where no imaginary 
part was find to that order. 

Some general arguments can be advanced, however, supporting that a nonva- 
nishing result should be found to higher loop order. Namely, in the second paper 
of [2] it was pointed out that (if it were applicable) the optical theorem relates the 
(technically quite difficult) computation of the imaginary part of the free energy to a 
much simpler tree level calculation, id est, the vacuum decay into identical particles 
[3]. A related phenomenon is the decay of a particle into several identical particles 
[4]. Besides, there is no reason for this effect to be restricted to de Sitter space; on 
the contrary, it would be natural to expect it to be quite generic. 

It is true, however, that all our intuition is based upon flat space examples, 
with the ensuing asymptotic regions, and S-matrix elements that can be computed 
through LSZ techniques. Outside this framework it is not even known how to define 
a particle to be decayed into nor the interacting vacuum |vac) in the absence of a 
well-defined energetic argument. 

A related issue is the study of the time dependence of transition amplitudes. The 
linear dependence in time is one of the key aspects of Fermi's golden rule. The fact 
that is problematic in curved space, where there is no naturally preferred coordinate 
system in general, has been remarked in [4]. It is to be stressed that use of non- 
cartesian coordinates is not without problems even in flat space-time, and this is 
even more true about polar coordinates in field space. One of the purposes of the 
present work is to examine this problem, by computing overlaps between states that 
differ by a finite time in whatever coordinate system we are using. 

The paper is organized as follows. In section 2 we review some results regarding 
transition amplitudes in both fiat and curved space-time, and we present the basic 
quantum-mechanical formalism that we will use in our calculations. It is based in 
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the functional Schrodinger picture for finite time intervals. In section 3, we put our 
techniques to work in order to recover the standard quantum evolution in Minkowski 
space-time. This we do in order to check our formalism, and to compare with or- 
dinary quantum mechanics (as opposed to field theory). In section 4 we apply this 
formalism to de Sitter (and briefly to anti de Sitter) to examine the (conformal time) 
dependence of its transition amplitudes. Some technical details have been relegated 
to the appendixes 
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2. Overlaps and evolution 

In flat n-diinensional space-time [5] (where energy conservation holds) differential 
transition amplitudes ^ from an initial state with Ni initial particles do behave for 
T large enough as 

where T is the time span during which the interaction is turned on, and Vn-i is the 



^It is not without interest to recall how the linear dependence in time appears on quantum 
mechanical survival amplitudes using old fashioned time-dependent perturbation theory (confer 
[6]). We start with eigenstates of a free hamiltonian 

Ho\k)=E°\k) 

and expand the full wavefunction as 

ip{t, f ) = ^ Ck{t)uk{x)e 



-iE"t 



Schrodingcr's equation 



then demands that 



which implies 



^l\4,) ^ (Ho + H,m 



-^eU 



ETT -iE°t ■ Y^ • -iE° 



Cfc = -*^Cfc,(fc|Hi(t)|fc>*(^'=-^'=')* 



k' 

When the initial state is an energy eigenstate of the unperturbed hamiltonian, id est, such that 
Cfc(O) = Skp; then, 

Ja 
where Ukk' ^ Ek — E^i . Assuming a constant perturbation leads to 



2 .i^^kpT , sin- 



c,{T) = ^ e^^ Hi 



Jkp 



T 



and the decay probability to an state |fc) is 

■ 2 ^kpT 

P,iT) ^ \c,iT)\' = 4 \Hl/ '-^^ 

In the limit when T — > cx3 this reduces to 

Pk{T)^\c,{T)\^ = 4\Hl/TS{u;,p) 

The survival probability is 

P(T) = 1- I p{k)dEkPkiT) 

where p{k) is the density of final states. 
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volume of the codimension one spatial sections of constant time of the system (so that 
in the particular case of vacuum decay it is proportional to the full n-dimensional 
volume). 

In the opposite limit, Maiani and Testa [7] have shown that there is a divergence 
at small times T — t- which survives even after renormalization. In order to elim- 
inate it, a careful study of the incoming wave packet is necessary. They studied in 
particular the example of an unstable scalar particle of mass M decaying into two 
other scalar particles with masses mi and m2. In the narrow packet approximation 
for the initial state they were able to prove that 

S{t) = 2'Kg^ e-**^* (ao + ai + /3i (t)) 

where M is the renormalized mass of the resonant state. 

The assumption will be made in section 2.1 that the quantum mechanical formula 

remains valid in curved backgrounds, where the hamiltonian is generically time de- 
pendent. 

2.1 Survival amplitudes. 

Let us introduce the general formalism first in flat space language, but in such a 
way that it is easily amenable to generalizations to curved space. The whole aim of 
the present work is to compute the overlap between an state |in) defined at a given 
time time tj and another state |out) defined at a different time tf (both times can 
be finite). This would become S-matrix elements in case tj — )■ — oo, t/ — ?■ oo and the 
interaction (including the one resulting from the background gravitational field, if 
any) is assumed to be switched off at asymptotic times. 
First of all, the survival amplitude is an overlap 

Mtf.u) = mf)\,p{t::)) 

If the spate is normalized so that 

then the unitary evolution (this is a crucial hypothesis) does preserve the norm, so 
that 

mf)mf)) = i 

Then Cauchy-Schwarz's inequality guarantees that 

\A{tj,u) I ^ \mf)mu))\ < \{i^{u)mu))\.\ms)mj))\ = i 
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This means that the quantity [T = tf ~ ti) 

T{T) = -^\og\A{tf,U)\ 

must be positive, and in case it is independent of T in the asymptotic regime, could 
be rightfully interpreted as the decay width of the state. We shall refer to it loosely 
as decay width even when it is not constant. 

Survival amplitudes are therefore powerful tools to detect instabilities; they are 
however somewhat blind to the final state of such decays; we will have no precise 
information on the decay products. This appears to be an important open problem 
from this viewpoint. 

The first principles path integral formula reads 

Sitf,U)\j= {ont\m)\j = J[D^f][Dcp,] ^t^[iffY{ipftf\ip.,U)\j ^t.N (2.1) 

where the integration measure [Dip] is defined in the space of field configurations at 
fixed "time". The wavefunctionals, which are functionals of this fields, are given by 

<i't^[ipi] = {ipai\m) (2.2) 

An external source J is introduced as usual in order to treat interactions by functional 
differentiation. 

The problem is then reduced to first computing the wavefunctionals of both states 
(itself a nontrivial task), and then the field transition amplitude, which is really the 
Feynman Kernel, or in modern parlance essentially the Schfodinger functional 

K[J]['^ftf,'^iti] = {yDftf\:^iti)\j (2.3) 

followed by a final functional integration over the possible values of the fields. 

This Schrodinger functional will be computed using the general expression in 
terms of path integral: 

K[J][^'t',y^t] = {y^'t'\^t)= f ' " ^^^e*^''^""^^^'^'^) , (2.4) 

J^{-,t)=ip 

where E^ is a codimension one hypersurface of constant time. 

Let us remind some well known facts in quantum mechanics (confer, for exam- 
ple, [6]). If the initial state |in, tj) (where we have explicitly indicated the possible 
time dependence) be it the vacuum state or otherwise, is an eigenstante of the full 
hamiltonian, then Schrodinger's equation imply that the modulus of the survival rate 
is equal to one and this is true for any values of t, and t/: 

\Ain{tf,ti)\ = |(int/|inti)| = 1 
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It is plain that the survival rate is nothing else that the self-overlap (|in) = |out)) at 
finite time interval. 

This means that the only way an state (vacuum or otherwise) can be unstable 
is by it being a superposition of energy eigenstates. Then the study of the survival 
amplitude for finite time is quite useful, because we do not need to know any details 
of the decay process (which is a complicated thing in the absence of asymptotically 
fiat regions). 

In order for a given state to be unstable it is not enough that the survival rate 
depends on time (this happens already for a linear superposition of only two energy 
eigenstates), but that this dependence has to be monotonic in time. It is enough, for 
example, that 

The actual dependence of the survival rate in quantum field theory with the time 
interval is however quite complicated. Besides the divergence at small times uncov- 
ered by Maiani and Testa [7] (whose understanding demands a careful treatment of 
wavepackets in the initial state), it can be explicitly shown in some models that the 
behavior is oscillating, except at asymptotic times (T = t/ — tj — )■ oo). 

3. Survival amplitudes in flat space 

It has been already advertised that the use of our techniques is best illustrated in 
the simplest fiat space example. It is going to be a rather long computation, so let 
us now draw its roadmap. There are three steps. The most important one is the 
computation of Schrodinger functional or Feynman Kernel; which is the quantum 
mechanical transition amplitude between states with well-defined values for the fields. 
This involves the computation of a determinant with Dirichlet boundary conditions. 
The final step is to integrate over the boundary values of the fields, weighted by the 
wavefunctional of the state, which we also need to know at this stage. 

3.1 Wavefunctionals 

This means that the first thing we have to do is to find the wavefunctionals [8] [9] of 
the states |in) and |out). Let us begin with the vacuum. First of all, in fiat space 
the free field and momentum operators read^ 

(x)= / -^^^^(a.^e'P^' + aU-'^^'^ 



nix) = -z I -J^^l^ (a,e^^' - a^^^^') (3.1) 



2 



^We will work in n dimensions, so the vector notation x, p, etc. will mean always an integration 
over 71 — 1 variables. 
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so that 

dx ^^ f 1 



in such a way that the vacuum wavefunctional obeys 



Up 



"^^ -e''''iv^pVix) + J-A^,\{v>\0) = 



(27r)V y ]l ujp6ip{x) 

so that the solution looks exactly the same as in [8] [9] , namely, 

^oM = (v^|0) = iVe-i^'^^'^^^(^'^)^(^)^(^) (3.2) 

with 

u (x, y)=u (f , x) = J T^^^e''^""'^'^^ 
Using the functional Schrodinger's equation 

i — ^o[v,t] = - dx i-^ - ri'Wiipdjip + m^ip^j ^ob,t] 

it is possible to determine 

with the vacuum energy defined as 

where the spatial volume is denoted by Vn-i = J dx. From the expression above it 
is plain that Eq is both ultraviolet and infrared divergent. 
The time-independent normalization factor is given by 

uj\l '^yi r — ^^ log ^ 



It is useful to consider eigenfunctions of the kernel defined such that 



In fiat space those are just plane waves 
and the eigenvalue reads 



e = ujp 
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In momentum space the vacuum wavefunctional reads 

where 

Then the one-particle state (defined in the non-interacting Fock space) would be 
defined as 

^i[0] = J dxUx){<P\l) = J dxUx){(j)\a\x)\Oh) = 

—^ / dxdyf,{x) i ^/uj{x,y) - ^"^/^(f, f)^Tpr 1 ^o[0] = 

VTe f dxMx)(Pix)^o[<P] (3.3) 



so that 



El = Eq + e = Eq + uip 



In the general case in which there are no asymptotically fiat regions of spacetime 
the uselfulness of those is quite hmited. 

3.2 Inclusion of the interaction 

First we write as usual the interacting kernel in terms of the free one using sources 

K[J] =e*4''^*'^^^^(il7)A'o[</] 
where Kq is the kernel that corresponds to the free action 

So[<P] -J<P^ j d^x (^-{d<pf - ^02 _ j(^)^(^; 

dt{-\4>k\'-^\<Pk\'-J-k{t)Mt)] (3.4) 



(27r)"-i4 V2'" ' 2 

In order to perform the functional integration, we follow Sakita [10] and split the 
field into a classical piece, (pKt) (with boundary conditions yet to be specified) and 
a quantum part, Xkif) 

0fc(t) = 0^.(t)+Xfe(t) 

so that we have for the measure Vcf) = Vx, as well as 

/dk /"*■' / \ 

j^^^ j^ dt [ci>lx-k - ujUlx-k - Jk{t)x^k{t)) 

(3.5) 
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The last piece can be written as 



dt ( ^ [4>lX-kj - X-k(pl - ^l4>lx^k - 3k{t)X-k{t) 



= <Plitf)X-k{tf) - <Pl{U)x-k{ti) - / dtx-k[(t^l + uji<Pl+3k{t)'i (3.6) 

and choose the classical field 0^ as the solution of the equation (pl + cu^^^ + jkif) = 0, 
so 



So[<P] -J<P = Som - W + Solx] + / rff [x(x)0^(x)] 



u 



(3.7) 



There is an additional contribution coming from the wavefunctionals in the sur- 
vival amplitude (2.1). For Fock vacuum wavefunctionals like (3.2), the exponent, 
depending on the boundary values of </>, can be written as: 



dk 



-^ujk (Iv^/fcP + \Vik\^) 



(3.8) 



2 J (27r)' 

I f dk / 

= -^J -^:;^^k[mU)\' + \Xk{U)\^ + 2x-k{U)(Pl{U)+ {tf term) 

The full monty of boundary terms in the sum of (3.7) and (3.8) is then 

i<Pkih)x~k{tf) - i4>l{U)x-k{U) - Uk {<pl{tf)x-k{tf) + X-k{U)<Pl{ti)) 

where the i comes from the one in front of the action, e*"^". This means that if we 
impose on the classical solution cff the boundary conditions 



i4>liti) + UkMU) = 



(3.9) 



this boundary terms vanish, and the classical field can be expressed in terms of the 
finite time Feynman propagator, At{x,x'), so that 



<P'{x) 



d''x'AT{x,x')J{x') 



Taking into account that 



^n 



dtdx 



dt 



d''x-{n + m^)(j)^-J(f)^ 



dx 



(3.10) 



j:^^mm' + \m.)f-)-'-f 
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(where the boundary conditions obtained in equation (3.9) have been used), the 
first terms cancel precisely with the remaining 10^^^ contribution in (3.8). The full 
classical piece in the exponent is given by an expression quadratic^ in J: 

J(j)^= r drxJ{x)(t)\x) = - r (f'x r rx'Jix)J{x')ATix,x') = -JAtJ 

Jti J ti J ti 

and we are left with the following expression for the free survival amplitude: 

AoitjMj = e^^' '^"^'^"^'^(^)^-(^'^')^(^'Uo(t/,t,,)|j=o (3.11) 

We still have to compute Ko[0] insofar as it depends on the initial and final times, 
ti and tf as well as the boundary conditions [ipt, iff], and integrate it convoluted with 
the remaining terms of the wavefunctionals. The exponents of these terms depend 
only on 

so the integration variables can be shifted: 



[Dipi][D^f] = [Dxl][Dxl^] (3.12) 



3.3 Classical solutions 



The way the computation has been organized is such that it stems from the careful 
evaluation of classical solutions with well-defined values for the fields at initial and 
final times. Let us first examine carefully the situation on flat space, and then build 
upon that. 

Any solution of the free Klein Gordon equation in flat space can be written as 

-—e'''''g{k)6{k'' -m^)= — — g-^'^"^-- {g+ cos ut + ig. sin ut) 

where 

_ 9{(^,k) ±g{-u,k) 
9± = ^ 

It is not difficult to show that the solution that reduces to 0i(x) aX t = ti and </'/(x) 
at t = t / = tj + T is given by 

_ dk _j^^ sinujk{t — ti)(f)f{k) — sinuk{t — tf)(j)i{k) 

X I — / 77: r- 7 C 



(27r)"~^ sinUkT 

The derivatives at the boundary are fixed and given by 



^The time integration in the definition of cp'^ takes place in the whole real line, while the time 
integration in the J(j)'^ term is constrained within the interval [ti,tf]. Nevertheless, the result can 
be proved to be independent of the value of the source J outside this interval. 



-11 



4>c{ti,x) 



dk 



(27r)"-i 



-ikx (pf{k)-cosoJkT(pi{k) 
e ujk — 



sinukT 



Mtf, 



X] 



dk „^^^ cosuJkT(t)f{k) - (f)i{k) 
e Wfc — 



(27r)"-i 



sin oJkT 



We will eventually be interested in the limit when T — )> oo. Choosing tj = — |- 



and tf = ^ it reads 



x 






-ikx 



sin cjfct cos cufct \ , .a / siii ^^fc^ cos cu^t 



sm -|- cos 2 



sm -|- cos -|- 



.(A:) 






-ikx 



""^)M-^.(*) 



sm 



'"i^ M+^.w 



cos 



(3.13) 



This formula does not hold when 
for a solution to exist is that 



0. In this case a necessary condition 



TT, 



ti-tf e -z 



u 



but this cannot hold true for all frequencies Uk- The most we can do is to make the 
solution vanish at one point. In this case, the field reads 



(/)c{x) 



d^k 

8^ 



36 '''^g{k) sin ujk{t-ti 



Let us examine the classical action when g = X = 



Sr 



d^'x { \ {d,<P^f - ^0^ 



1 



1 



rx-d^ (0e5^0c) - 7,<Pc (n^c - ^202) = / dx- 0,0, 



dk Uk 



(27r)"-i2sinwfcT 
In the particular case when T — )■ ioo we get 



{i\i<Pf)k\^ + \i<Pi)k\^) coswfcT - 2Re[(0j)fc(0,)_fe]}(3.14) 



Sr 



-i I dx- drda^ ' 



'^c^OVd 



i dkujk {(f) f {-k) (j)f (k) + (pi {-k) 0j (k)) 



3.4 Computation of the Schrodinger functional. 

Let us compute the Schrodinger functional i^o[0][Xi5 X/] that is, the transition am- 
plitude between states with well-defined values for the quantum fields in the free 
case. This is the quantity that in norelativistic quantum mechanics is aptly named 

Feynman's kernel. 



Komx^tf,x'U]^ I Vxe 



i5oM 



(3.15) 
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This functional integral is computed with fixed Dirichlet boundary conditions at the 
endpoints and its functional form is: 

Ko[0][x^tf,x'U] = e'^^['^''^^ldet(n + m2)-i (3.16) 

This expression can be achieved expanding again around a classical solution, 
given by: 

X(t,f) = / , ^,^, ^^ ^^"^^^^^ " ^^ + ^^ ^^"^^(^ - '^ e-^'^' + at, X) (3.17) 



(27r)"-i sinwfcT 

so the classical action read 



1 r dk 



f 



S=lx-, X'l = 2 y p;^ Xrt-. . = (3.18) 



and the only thing that remains is to compute the determinant 



det(n + m2)=/ r'^e'^ot^l (3.19) 

Jo 



with vanishing boundary conditions U, = ^\, =0. 

It is quite easy to check that the eigenfunctions are given by 



^^ = V T ['''' T ) ' 

with eigenvalues \j = ul — ^r^, where j = 1, 2, . . . For k = they are normalized in 
such a way that 

/ dt^n{t)^m{t) = Snm 

Jo 
The associated zeta function (confer [11] whose coventions we follow here) reads 

^.^_^^ f dk ( e + m^-^^ '-' 
^^'^ = ^"-^ J M^ ^ [ 7 

where Ki-i = / dx. Now we can use the following identities 



(27r)"-i ^ ^ (27r)"-i r(s) 



oo 



X-'T(s) = / dr t'^' e 



-1 „-At 
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so that 






n^^^^V^.^i-ir^^ ^ r ,^^s-^i ^-v - 






2sT/ ('^2\-s+iV^ 1 /-oo 



(27r)"-ir(s) 2 io 



6*3 ( 0, 



(3.20) 



where we have performed analytic continuation in the mass and Q-^i^z^ r) is the ordi- 
nary Jacobi elhptic theta function, 



oo 



■^^ „2 



-Kmz A,'n?"KT 



d^{z,r)^ ^ e^-^"V 



n=— oo 



This expression is divergent at s = 0, precisely the point at which its derivative 
is wanted. We can remedy using a Poisson resummation in the first summand, id 

est, 

^ / itTT \ mT ^ ( im'^T^ 
9s 0, — - = ^=^3 0, 



and representing by l? = — ^"^^2 , 



n— 1 . o , n — 1 



«^'-'») - ^IS^"^'-"-' {y^f--'-'-^' (°"^) - 2^^ -) } 



The integral is now convergent for appropiate values of s: 



62t; Vl?!/ ^^ 2-- V 5 y ' V 2 

bl>i 



and then 



n — 1 . r). n — 1 



C(5) 



2(27r)"-ir(s) 



/i"K-i |^(!!I^^tZi^ E^'^-^^t-^ (2^-^^') + (3-22) 



imT / n\ ( \ — n 



For even values of n: 



n-l 



^-^^ ' ^^^1^1^ [7 - log(mV/x^) + ^(°Hl + n/2)] ^ (3.23) 
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(3.24) 



The sourceless Feynman kernel is precisely 

^o[0][x^,xi=e^^=e-K'(o) 

3.5 Flat space survival amplitude 

Let us study the vacuum survival amplitude in the free case performing the relevant 
integrals by brute force, id est, without introducing any sources for the fields 

Ao{tf,U) = {Otf\OU) = f[Dipf][Dipi]{Otf\ipftf){ipftf\ipiU){ipiU\Oti) = 



where we have introduced the vacuum wavefunctionals of section (3.1), and the free 
Feynman kernel has just been shown to be 

Ko[0][^ftf,^,t,] = {iPftf\ip,U) = det (D +m2)-^ e'S.i^.,^,] 

where the determinant does not depend upon the boundary values for the field vari- 
ables. 



This yields 

5o(t/,t.) = liVpdet (D + m'Y^- I [Dipf][Dip, 



■exp 



dk 



(27r)n^i ^^ ('^A ^*k) 



\ cot oJkT 



2sincjfeT 2 ' 2 

We shall dub the functional determinant of the operator 



2 sin cjfeT 

+ \ cot uokT 



B 



+ \ cot OJkT 



1 



2 sin tjfcT 
1 
2sintjfer 2 ' 2 



1 



^ cot ujkT 



the boundary determinant. 

The eigenvalues of this matrix are: 



A 



ztan(a;fcT/2) , 



zcot(a;fcT/2) 



Vfk 



(3.25) 



so the zeta function we have to consider in irder to compute the boundary determi- 
nant is (recovering the Uk factor we dropped in these eigenvalues): 



C(s) = K-i/i^ 



dk /Uk 
(27r)"-i \Y 



oi'^kT/2 



cosuJkT/2 



Vn-lH^^n-2'm 

2"^(27r)"- 



n—l—s poo 



dx(x^ 



1) 2 X 



1-s 



f.i^kT/2 

sin uJkT /2 

ixmT/2 



COS xmT /2 



(3.26) 



ixmT/2 

sinxmT/2 
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Now we can perform use the following expansion: 

where 7ro(s) = 1 and tijj^i{s) = s ■ . . . ■ [s — j) . The zeta function is then: 



O— s|'2yT-yt— 1 I I \ 1 _|_ p—imTx j I 1 p—imTx 



-2 



E^^^.(^) (3-28) 



(27r)"-i ^ 2j! 

with /2j(s) = ^^ dx{x^ - l)^xi-^e-2^™^"^\ 

Since we know also that 7rQ(0) = 1 and 7r'^]^(0) = (— l)-'j!, we have: 



l^n„_2m'-2 I /;(0) - log(-m/,,) /„(0) - J^ Yj'^AO) \ (329) 



The integrals can be calculated then: 



/2,(0) = -^r (^) (2jmT)i-tir|(2^jmT) 
In the even n case the zeta function corresponding to the boundary determinant is 

(3.30) 
Collecting the results of this paragraph with the ones of the previous one, the 
vacuum survival amplitude in the even n case reads 



Aoitf.ti) =exp 



1 7r 2 m" ^ { f 1 — n 



2 2(27r)"-i t V 2 

(3.3i; 



^ (.mTKlj^ _ i^g(^2/2^2) ^ ^(o)(i + ^/2)] 
yTi\n/z,y. 



In is remarkable that the end product of this computation is of the form 



v4.o(t/, tj 



oc e 



-iEoT 



16- 



with the vacuum energy given by 



En 



1 r 

1 TT 2 rn 



;-i)^ 



K-i 2 2(27r)"-iv^(n/2)! 



[7 - log(mV2/i2) + ^(0) (1 + n/2)\ (3.32) 



The inclusion of the interaction in these considerations can be achieved through 
Feynman diagrams with finite time propagators built in them. 

3.6 Adding particles 

If we try to perform the same calculation above for an excited state (id est, an state 
containing particles), the only change is the boundary wavefunctions in the second 
determinant: 



[D^i][D^f]>^i{py^f{p) exp 



dk Uk , . f i cot ujkT — 1 —i esc ujkT 

(2^0^^^ ■^ V -i CSC uJkT icotuJkT-1 

(3.33) 



Vfy k 



for the case of an state with a particle of momenta p. 

If we introduce sources coupled to the boundary values of 0, the calculation only 
needs the addition of a term coming from the derivatives: 



[Dipi] [Dipf] exp 



dk 



(2vr) 



- {^^ Vf)^kMk 






+ {J\J^)^ 






(3.34) 



exp 



dk 



(27r)"-i^ '-'' *= 



jf 



[Difi] [Dipf] exp 



dk \ 1,^ I '^i 



We can get then the additional factors due to the presence of particles: 



k 



ujk \e 



6 



SJ'U6Jl 



exp 



dk 



4 



4y(2^(''^-4)A'r 4 



oc e 



-icOkT 



(3.35) 



(3.36) 



j=o 



and this means that the energy in the exponent e~*^°^ calculated for the vacuum- 
to- vacuum amplitude increases precisely in Uk- The reason why we did this check 
is that the way the linear dependence in time appears is quite different in quantum 
field theory and in quantum mechanics, as we detail in the next paragraph. 

3.7 Quantum mechanics 

It is useful to contrast the field theoretical calculation above with the quantum 
mechanical harmonic oscillator [12] 



ff 



-17 



The vacuum survival amplitude for the harmonic oscillator with unit massn is 
given by 



C(T) = (Ot = T|Ot = 0) = /"l)g,I)g/e-'^«'/2g-a;g2/2 f ^ jyq^iS[q] 

J Jqi 

Expanding the trajectory around the classical solution, q = qc + y, we have: 
C{T)= /"pg.Pg^e-'^'/^e— 9?/2giSfe] f Vy e'^^^^ 

J Jy{0)=y{T)=0 

The action for the classical trajectory can be expressed as: 



(3.37) 



(3.38) 



Sh 



1 



f 



u 



2sinco'T 



{qf{cosujTqf - qi) - qi{qf - cosuTqi)) (3.39) 



and this is quadratic in the boundary values for qc so that the boundary integral is 
gaussian: 



/"pgiI)g/e-"'?'/2g-..9?/2gi5fe]^ ^ 



- 1 + I w cot uT — I w CSC uT 
— |w CSC uT ~ f + l^^ cot uT 



(3.40) 



-1/2 



71 { -coUl-i 



V2 



J CO 



The Dirichlet determinant can be computed through discretization, following Feyn- 
man's original argument, to be 



det(--+c.^ 



CO 



27!-ismuT 



1/2 



It is amusing to remark that a zeta function computation gives this same determinant 
only up to a constant. Altogether it yields 



C{T) 



e-f-^ 



The differences with the quantum field theoretic computation are now clear. The 
dominant terms for large T come from the Schrodinger functional only in the field 
theoretic case, whereas as a result of cancellations, they come from the boundary 
determinant in quantum mechanics. We just checked that nevertheless the energies 
of the excited states are also correctly given in quantum field theory. 
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4. Survival amplitudes in de Sitter space 

Let us now turn to the main object of our interest, namely the (in)stabihty of the 
vacuum state in (anti) de Sitter space. We shall mainly use here the de Sitter metric 
in horospheric (Poincare) coordinates, where z plays the role of conformal time 

as = :^ 

z^ 

The conformal time is positive semidefinite 

0<z<oo 

It is sometimes useful to write z = e^^*, where the Hubble constant, H is related to 
the radius by if = y, so that the metric appears in the steady state form 

ds'^ = dt^ - e'^^^5ijdx'dx^ 

In these coordinates it is plain that in the limit if — )■ (/ — )■ oo) flat space is 
recovered. 

We are interested in the survival amplitude of a certain state |in) between (con- 
formal) time z and z' (both times can be finite) 

Ainizf,Zi) = (inzflmzi) = / [L'v?/] [i^v^j] {in z f\ip f Zf){ipf Zf\ipiZi){ipiZi\mZi) 

4.1 Wavefunctionals 

We shall expand the free field as 

(f){z,x) = / dp{apVp{z)e'^'^ + alvl{z)e''^'^) 

and the canonically conjugated momentum 

dp 



TT[Z,XI 



{apvlXz)e'^^' + al{v'^nz)e-^^^^^) 



Iz'^-^ 
Our modes are normalized by the usual Klein-Gordon invariant scalar product 






The creation and annihilation operators are given by: 

ak = ~i—j— / dxe"'^"" {{v'f,)*{x)(j){z,x) - Iz''~^vI{z)tc{z,x)) 



t -^ 
al = i- 



2-n 



I 



^^g+ife [y'^^x)(t){z, x) - Iz^'-^kizMz, x)) (4.1) 
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It is now quite plain (at least formally) how to compute wavefunctions for dif- 
ferent states. Let us begin with the wavefunction of the free (Fock) vacuum. It is 
defined for appropiate destruction operators and a given conformal time z by 

OfclO) = I dxe-'^^ {y'^{zY<\){z,x)-vu{zflz^-'^'K{z,x)) |0) = 

we are thus led to a differential equation common for the vacuum wavefunction, 

dxe-'^^ {iVk{zYz''-H-^^\v'^{zrif{x)\ (v9^|0) = 

The vacuum wavefunctional (92 10) is the exact analogue of the Schrodinger wavefunc- 
tion ip{q,t) = {q\ip), where the completeness relationship ^ \q){q\ = 1 is assumed in 
a time-independent way. Here we introduce a time-independent basis \ip) such that 



This basis is defined in such a way that the diagonalize the field operator 

at a certain fiducial time, z. But the basis itself depends on this fiducial time in 
a nontrivial way, and this we have attempted to represent by writing explicitly the 
basis as {(f)z\. It follows that a gaussian ansatz 

(mz\o) = ]\f Q-h I '^^'^y^^(^'yM^)v>(y) 
is indeed a solution, provided 



(27r)"-i J Iz''-^ Vp{z)* 

This gives a natural definition of non-interacting vacuum state corresponding 
to the modes Vp{z). The present definition of vacuum depends on the modes used, 
and this in turn depends on the physical setup of the question asked. This is a gen- 
eral problem of quantum field theory in a curved space, not specific to our formalism. 

Through the functional Schrodinger's equation perturbative corrections to the 
noninteracting vacuum can easily be found. The concept of particle is a delicate one 
when asymptotically fiat regions are absent. A possible definition of a multiparticle 
state in the present context is, for example, 

{<^ z\ki . . . kp) = {^ z\al^ . . . aljO) 

but it is plain that the usefulness of such a definition is quite limited. 
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4.2 Classical solutions 

The action for a scalar field in a generic conformally fiat space (of which both de 
Sitter and anti de Sitter are particular instances) can be written in a very simple form. 
We shall insist for no particular reason in keeping the coordinate z dimensionless, 
so that the dimensionful coordinates are x^ = {x^,x^, . . .x"'~'^,lz). The metric is 
conformally fiat 

ds'^ = a{zY rj^i,dx'^dx^ 
In de Sitter space the coordinate z is timelike, so that it follows that 

S,s[^] = I Idz dx a" (^1 {^,{dAf - (V0)^) - ^0^ - ^0^ 

In anti de Sitter, owing to the fact that the z coordinate is spacelike, this reads 

SA,sm = Jldzd-^'x a" Q-^ (-^(5.0)^ + 0^ - (V0)2^ - ^0^ - ^0^ 

(where now x includes all coordinates except x^ = t and z). 

We shall actually redefine the quantum field (but keep the same notation for it) 
in order to shift all depence on the background towards the potential 



^new = a 2 0old 



The lagrangian now reads (remember, now = 0ncw) 



c^lio,y-'^,'^Y^^-^^,'± ^,, ^ 



where 



,2/\ _ _2„2 



n\ a (n \ (n \ a 



2 



m'(z) = m^a' ±1--— r^--2--l, ,„ 

A(^) = a^-" A (4.2) 

(where d = ^; the upper signs are for de Sitter space, and the lower ones for anti de 
Sitter). 

In both de Sitter and anti de Sitter, a = - which gives, paying due attention to 
the fact that z is dimensionless, 

S,s = \ [dxldzU {d^<pf-{V<Pf- '^'^'~ 7^ 0^ - ^'^A+ fdx''-^-' 



2 J 1/2^'^' ' ^' P z^ ^ 12 I ^ \ J Azl 

Sa,s = I / d-'xldz -i (9.0)^ + 0^ - (V0)^ + /^, ' 0^ - ^^0^ 



^/ 



21 



where care has been taken to keep in de Sitter all boundary terms for future use. 
Incidentally, those are totally irrelevant for anti de Sitter, because we are only inte- 
grating the time variable over a finite time interval; but they are quite important for 
de Sitter space, because they enforce a change in Feynman's propagator as explained 
in the appendix in some detail. 



It is amusing to remark that up to a constant factor the 7 factor defined in the 
3ndix as 7 = f^^ is just the de Sitter temperature T = ^ 



n-2 ^ 

7 ?T2 = 27rT 

z 

It is well known that this temperature is associated to the unavoidable presence of a 
horizon because of the lack of a globally timelike Killing vector [13]. This coincidence 
is due to the fact that there is a single energy scale in de Sitter space. 

This action can (and will) be interpreted as a Minkowskian action for a massive 
field, with a time-dependent potential given by 

The perturbation is a time dependent one for de Sitter space; whereas space de- 
pendent in anti de Sitter. In this split between free and interacting hamiltonian all 
information on the curvature of the space has been dumped into the potential term. 

4.2.1 A different split between free and interaction terms. 

It can be more convenient for some purposes to keep all quadratic terms (dimen- 
sion two operators) in the free lagrangian and treat higher dimensional operators as 
interacting terms. 

The scalar field action in de Sitter space 



'■^ Idzdx ( z^ 



reads in momentum space 

/dk 



m 



^f Idzdk f z"^ ,.„ . ,o ,o.o, , .n\ m^ 
Writing as in flat space 



+ Xk 
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the equation of motion reads 

9k 9k + 2 9k = 

whose general solution reads 



To find the solution that reduces to 0i(x) at z = zi and to 02(5') at z = ^2, let 
us define 



</l,2 = '^(^\^^m?-p. iklZi^2) 



Then 





n,2 = V(„-i)2 


mn^ (^^^1,2) 


Ci 


1^201 (A;)^i=^ - 


-Fi02(A;)^2' 


^1^2- 


-J2YI 


C2 


-J2Mk)zi^ 


+ JiMk)z2' 


J1Y2 


-J2YI 



(4.3) 



(4.4) 



and 



0f -(2) = Ci J, (A;/z) + C2Y, {klz) 
where u = (^^) ~ "m^/^. For any classical solution the action on shell is given by 



Sc = I dx-^<pd-A\i = (27r)'^-^ / - -—^^,(z){d,<Pu){z 



1 



Let us now find the solution of the equation of motion with a delta-function 
source and Dirichlet boundary conditions. 
Using the fact that the Wronskian 

W[Uz),Y,{z)]=2 



ZCOSUTT 

the Dirichlet propagator is easily found to be 

Af (z, zo) = ^ e-^"" klzo9 {z - z^) costtz/ {-Y, (klzo) J. (klz) + J, {klzo) Y, {klz))+<pf^' {z) 

We shall develop the expansion associated to this propagator in a forthcoming 
work. 
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4.3 The Schrodinger functional 

The Schrodinger functional is given by finite-time Feynman's diagrams with position- 
dependent vertices 



A{Zf,Zi)\j = e"'-^^'/''^^'^^^(^'''77) g|/z7<^"^'i"^'''W^T(2^,:r')J(x') 



-^0 1 J=0 
J=0 



To a given order in perturbation theory, it corresponds to vacuum diagrams (the 
same that contribute to the usual vacuum energy) computed with finite time Feyn- 
man propagators. The important thing to notice is that the only dependence on the 
boundary values of the fields stems from the classical action. 

The first diagram to be computed is the "circle" , which is simply: 

r^ f dk 

Mo,o = / Idz dx-——^/\T{k)[z,z] (4.5) 

It is to be remarked that even this diagram carries some information about the cur- 
vature of the space through the 7 terms in the propagator. 

We can take advantage of the specific form of the said propagator, in the sense 
that only the first coefficient (C.13) contributes to the simple diagrams we will con- 
sider. In this case, the amplitude reads: 

(4.6) 
where we have neglected the product 7j7f. 

In the limits of large and small Z, (physically, the relevant quantity is Zml) we 
have: 

Mo,o^°^^1^^{-<^.o).|-,z=.(^,-:)} (4.) 

where 

/oo 
dxix"^ - ly x'' 

roo / -2imlZx _ i \ 

MZ)^l <i^-(^-=-l)-L__^ (4.8) 
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The second diagram is the same as before, but with a mass insertion: 



Mi,o 



Idz / dx- " " ^ '^ ^ ' 



'j^:;^i^Tik)[z,z]ia + 



(4.9) 



with q; = —771^/2 and /3 = (m^ — n{n — 2)/4/^)/2. We have a part proportional to 
the first diagram, and a second part proportional to: 



Zfl 



-ItOke-'^'^'^'f^''^ 



— 2ilmukZfZi {'yf(Ei{2ilzfUk) — Ei{2ilziUk)) + 



+-fie^'^'^''^'f+'^\Ei{-2ilzfUJk) - Ei{-2ilziUJk))) - me^'^'^'''^{-ffZf + -fiZi) + 
+ e^''''^'f{m{-ijZf + 7,z,) + 2iuJk{zf - z,)) 

where we have neglected again the terms quadratic in the 7's. 



(4.10) 



For large Z this contribution is just a constant independent of Z, while that for 
small Z has a linear and a quadratic part. In this last limit, the full amplitude reads: 



Mi,o 



Z-5>0 



a + 4 Moo + 



n-2 



ZJ 



2; ^^^00 



Vn^l^n~2m 

2(27r)""i^f 



ilZ^J 



n — 3 



(4.11) 



where the Mqo has to be understood as the small-Z limit shown above. 



The third diagram (the first contribution of the self interaction) is given by 

(4.12) 



Mi,o = 3t I Idz I dxj^r^£^AT{k)[z,z]AT{p)[z,z]X{z) 



(2vr; 



2n-2 



with X{z) = z"' ^A/24. The diagram is then proportional to: 



Idz — ■ — 



7im (-1 + e2^''^'=(^'-^') j + 2iuk) (lim (-1 + e^'^''p^''"'A + 2iUp) e^'K+^p)(2^'-^/-^0 



-2i/a;fcC^pe-^'('^'=+^^)(^-f+^>) 



7,/^m23-"a;^a;p(i/a;fe)-"-^e2^'(^/(^'=+'^'')+'^*^')(r(n - 3, 2^/^/0;^)- 



- T{n - 3, 2ilziC0k)) + -fil^m2''-''cUkto^{ilcOp)-''-'e^'^^'f^''>'+''-^+'''''^\T{n - 3, 2ilzfCUp)- 

- r(n - 3, 2Uz^Up)) - ^- [z]-"^ - z^-^) e2*'^/("'=+'^-)(^(7/ - l^){ook + oo^) + 2tUkU,)- 

It o 

- -ffm2^-''up{-ilukf-''e^'^'^''^f(T{n - 3, ~2ilzfUJk) - T{n - 3, -2ilziUk))- 

- 7/m23~"wfc(-i/a;p)^-"e2^''^'=^-f (r(n - 3, -2ilzfUJp) - T{n - 3, -2ilziUjp)) (4.13) 
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Mo,i 
Mo,i 



For large Z, the leading term is proportional to Z" ^ 



,2 _2n-4 f /^_3\2 ^^n-3 



32(27r)2«-2 

Z^O '^^i^n-2 



iAfi2 v;-im2"-4 



32(27r; 



2n-2 



J 



/J(^,o 



n — 3 



+ 



^'^•°hS«^' 



/Z 



n-4 



+ 



-2, 



2 

ra— 5 



-^rv('^,o)z+ 



(4.14) 



^7*^ ( -^^, -1 - (n - 4) J ( ^^, 



n — 3 



n — 3 



^—2imlZx^ 



where: 

/'CxD /'OO 

hi{Z)= dx rfy(a;2-l)V(y2_i)V 

Ji Ji a;(2a;|/ + 7i(x + |/)) 

4.4 Vacuum Wavefunctionals 



(4.15) 



When using, as we do, the variables (pnew the appropiate starting point for the vacuum 
wavefunctional is the Minkowski one. 

The functional Schrodinger's equation, which stems from our main hypothesis on the 
Feynman kernel in curved space reads 

It is possible to solve it in a perturbative way in X{z) and A = m'^{z) —rn^ by writing 



^[0] = ^o[0] + A(^)^io[0] + A(^)^oi[0] 



In this way 



^ 



10 [ 



0^ 



(i^. 



1,22; / {K^^cp^y 



(4.16) 



as well as 

2 

It is remarkable that both terms are proportional to ^o; so that the total vacuum 
wavefunctional can be written as 



oi[<, 
— i 



yJ^vwf'u 



lYT[h-n,iz / {K, 



(4.17) 



^[0] =^0(1 + A5i^ + A(52^) 

It is also possible to view the functional Schrodinger equation as an evolution equa- 
tion, and assume that at a given conformal time z = z^ the wavefunctional is ^E'o[0], 
and then compute its future evolution in the conformal time. This is not what we 
have done here. 
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4.5 Survival amplitude 





\2 



m[z) 



Order (1,0) 

Order (0, 1) Order (1, 1) 

Figure 1: The first few diagrams that contribute to the vacuum energy. 




The only step still left in order to compute the (vacuum) survival amplitude 
is the integration over the boundary values of the fields, weighed by the vacuum 
wavefunctions as well as the classical action. We know already from our previous 
computation (confer equation (3.32)) that this contribution is subdominant in the 
large Z limit, and besides it preserves the modulus of the (exponentiated) survival 
amplitude, so that it gives vanishing contribution to the width . 

There are however calculable interaction dependent corrections to the vacuum 
wavefunction (as to any other wavefunction) ; they can be obtained through the 
functional Schrodinger's equation to any given order in perturbation theory along 
the lines of the subsection 4.4. We have not attempted to compute the effect of those 
corrections on the width. 

In conclusion, the value we get for the width of the vacuum state in the asymp- 
totic regime Z — > oo under the approximations of the present work is 



T{Z) 



^ 2(27r)"-iZ 



Re/oo(^) 



AK-inL2---S,^._,j^,^^(^^ (4.18) 



16(27r; 



2n-2 



5. Conclusions 



In this paper we have concentrated in computing overlaps between arbitrary states 
(in particular the vacuum) defined at two different times such that they span a finite 
time interval (were this interval infinite they would become S'-matrix elements, in 
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case those happen to be well defined). This has been done because there is some 
initial doubt as to how to define the good observables (id est, the analogous to the 
decay rate, for example [4] [2] ) which would presumably involve some sort of square 
of the overlap matrix elements themselves. 

The most important quantity we have analyzed is the survival rate, or self- 
overlap at finite (conformal) time span. This in turn determines a decay width in a 
straightforward way. It is found that there some effects already at tree level, which 
are presumably related to particle creation in the presence of an external non-static 
gravitational field, but we have not idetintified them unambiguosly. Our compu- 
tations are consistent with them being transients. They are however of potential 
physical relevance in the physics of the infiationary epoch. 

At the next order in perturbation theory, there is a new contribution which de- 
termines the vacuum width in a precise way in the adequate spacetime dimension 
(namely, n = 5 were the dependence of both Ioi{Z) as well as Ioo{Z) on their argu- 
ment subdominant). Further diagrams should be studied before a definite conclusion 
can be drawn on the main issue. 

This computation has been done for a particular wavefunction, which does receive 
corrections owing to the interaction. Other states can easily be studied within our 
framework. To the extent that fiat space computations are a good guide, we do not 
expect those improvements to change the physical picture dramatically. 

All the physical quantities studied in this paper turn out to be observer-dependent. 
It is not completely clear what could be the physical meaning of some phenomenon 
which is coordinate dependent (or what amounts to more or less the same thing, 
observer dependent). There are by now many examples of observer-dependent phe- 
nomena even in Minkowski spacetime (of which the Unruh radiation [14] observed 
by an accelerated oberver in the Minkowski vacuum is perhaps the best known); this 
does not necessarily mean that their physical meaning is fully understood. 

On the other hand, it is well known that the usual semiclassical approximation to 
the full quantum theory of the gravitational field interacting with arbitrary matter, 
namely quantum field theory in a external gravitational field treated classically (upon 
which the latter identification is based) is only an approximation to the true equations 
of motion, to wit 

(vac|^— ;[^^j„V^i]|vac) = 

Where the total action is the sum of the Einstein-Hilbert part depending on the 
metric only, the matter part, which depend on the matter fields, denoted here col- 
lectively hj ipi i = 1 . . . N, and the necessary counterterms, which depend on the 
metric as well as on the matter fields. Including sources, 

S[Jf^u, Ji] = Seh[9] + Sjnatt[g, -ipi] + Scount[g, V^i] + « / d''x^/\g\ I J^'^g^u + ^ ^Vi 
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The equations of motion are always formally true because they can be written in 
terms of the full partition function 



Z[J^i,,Ji] = / Vg^t^YyVipi e 



f'J ['-' i-ii^ I'-'i 



as 



5S 



Sg 



111/ 



6 1 6 



^ \'J iivi 'Ji_ 







j=o 



i 6J^u{xy i 6Ji{x)_ 

It remains to give a working definition of the composite operator g'^'^, but at the 
perturbative level this can be done. The state |vac) is the one obtained through the 
boundary conditions imposed on the path integral. 

The semiclassical framework states that this vacuum can be approximated by 
the matter vacuum in a fixed gravitational background g^i, 

I vac) ~ |Omatter)g 

This can be proven to be the dominant term the first term in a 1/A^ expansion [17] 
of a theory of gravity interacting with A^ identical matter species, but it is difficult 
to believe that this is the only instance in which this semiclassical approximation is 
physically reasonable. A general analysis of its validity would be welcome. 

Observables in the full quantum gravity theory should presumably be gauge 
invariant, that is, diffeomorphism invariant, and thus independent on the observer. 
What seems to be needed here is a gauge invariant definition of vacuum decay. 

More comprehensive computations are in progress taking into account the dy- 
namics of the gravity sector. 

Acknowledgments 

This work has been partially supported by the European Commission (HPRN-CT- 
200-00148) as well as by FPA2009-09017 (DGI del MCyT, Spain) and S2009ESP-1473 
(CA Madrid). R.V. is supported by a MEC grant, AP2006-01876. 



-29 



A. Flat space vacuum stability 

Let us first review the reason why those effects vanish in fiat Minkowski space, 
deriving in this way a simple formula for it that subsequently could be applied to the 
spaces of our interest in different physical situations. To assert that the fiat space 
vacuum state is stable is equivalent to assert that the free energy is formally real 
in Minkowski space, which through the optical theorem ensures the stability of fiat 
space versus multiparticle decay. Using LSZ reduction, the S-matrix amplitude for 
vacuum decay to four identical particles with wave functions u^^ . . . u^^^ reads 

(A;iA;2A;3A;4|0) = i / V | (7(1/1 )M"|/i<(l/i)V|^(l/2)M"l/2< (1/2)^1^(1/3) M"l/3<(l/3) 
Vl^(l/4)M"l/4<(l/4) [^y, + m^) [Uy^ + m") [Uy.^ + m^) (d,^ + m") x 

(O|T0(l/i)0(l/2)0(l/3)0(l/4)|O) 



At tree level 



A 



(O|T0(i/i)0(i/2)0(l/3)0(l/4)|O) = -j (ry^\M\^ {yi - y) A (1/2 - y) A (1/3 - y) A {y, - y) 
where the Feynman propagator obeys 



(n + m") A(a: - y) = —===5{x - y) 

\/\9{x)\ 



In fiat space 



—ikx 



Uk = — , ^e 



so that the net output at tree level is 



/ 



r^> , ^ _ p-^(fci+fc2+fc3+fc4);/ _ (27r)" , ^ — ^ (ki + k2 + ks + kA 

^UJlUJ2UJz^A ^UJiUJ2UJ^UJi 



This implies in particular a delta function on the sum of all energies. 



which does not enjoy support on physical particles. 

At the same time this gives a simple condition (assuming LSZ reduction is still 
valid) for this amplitude to be nonvanishing in an arbitrary spacetime, namely, the 
vacuum is unstable with respect to decay into four particles whenever the fourfold 

(9("=4)['0j overlap, where 

0(") [0] ^ f d^'y^\M\nk, (y) • • • Uk^ {y) 
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has got nonvanishing support on physical states. It is plain that this depends on the 
value of the determinant g (there is always a gauge in which (7 = 1) as well as on the 
set of modes Uk- Another quantity of interest in conection to a single particle decay 
into two or three identical particles is 



B. A first (naive) look at overlaps in de Sitter space. 

Let us examine the overlap with several different particles (id est, different coordinate 
systems) and in different spaces assuming LSZ reduction. The purpose of the present 
paper was precisely to improve upon this analysis, which we want to briefly present 
here. 

To begin with, let us assert that there is no vacuum decay, nor single particle 
decay into two or three identical particles in static coordinates (which exist for both 
de Sitter and anti de Sitter). We call static coordinates ones adapted to the timelike 
Killing, in such a way the timelike coordinate is ignorable. There may be many 
different such systems for a given spacetime. 

The reason is in them the exact modes have got a piece 

Uk ~ e-'^'fk 

where the functions fk do not contain the variable time. This is enough to produce 
a delta function 

which do not have support on positive energy particles. Using the fact that 

a;^ - m^ > 

it is also possible to show that there is no single particle decay in the static case. 
To be specific, de Sitter in static coordinates reads 

ds^ = (1 - ^) rft^ - ^ - r'dnu 

where the radius of the spacetime is related to Bubble's constant by 
and anti se Sitter space in the same coordinates, 

ds^ = \i + -^\de - ,^-rHnl_, 
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In anti de Sitter space the horospheric coordinate is spacelike, and it will be denoted 
by X, so that Poincare coordinates admit a FRW-like form 



ds^ = -dx^ + e^ (dt^ - ^ dyf\ 



which is manifestly static, so that there is no vacuum energy decay here. 

De Sitter space in global coordinates (spherical spatial spacelike sections) reads 

ds^ = dr'^ — coshT'^dQ'^_^ 

Again, in anti se Sitter in global coordinates, the metric reads 

ds'^ = cosh'^ rdO'^ — dr'^ — sinh^ TdVt^_2 

Finally, when hyperbolic spacelike sections are considered, de Sitter metric reads 

ds^ = dr"^ — sinh^ r {dip"^ + sinh"^ ipdVL'^_^ 

In contrast, anti de Sitter space in the same coordinates yields 



ds = sink xdtp — cosh x^X ~ cosh xd^ 



n—3 



which is again explicitly static. 

Lest the reader has the impression that anti de Sitter looks atatic in all coordi- 
nates systems, let us mention stereographic coordinates, in which the metric reads 

ds'^ = Q'^rj^^dx^dx'^ 

and for de Sitter space 

1 

n 






with x^ = rj^yX^x^ = t'^ — r'^. For anti de Sitter space 

1 



n 



i + S 



In this coordinates, the global staticity of anti de Sitter space is not manifest. 

First of all, a general observation [15] [16]. The link between the field modes and 
the particle concept is through the WKB approximation. Indeen, the Klein-Gordon 
equation 

yields, for 

6 = e^f +■■• 



the mass shell condition 



g'^'d^SdpS 



m^ 
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encoding the definition of positive and negative frequencies for the solutions at hand. 
By covariantly deriving the expression above we get 

V^d^Sg'^^dpS = 

Now, for any scalar, 

(V„V/3 - V/3V«) S = 

it follows that 

g'^^dpSV^V^S = 

namely the geodesic equation. This means that the vector -u" = ^^ is the tangent 
vector to a geodesic, which in turn implies that the hypersurfaces S = constant are 
geodesic orthogonal. 

C. Finite time propagators 

The general solution of the Klein Gordon equation 



[a + m^)(i) = J 



can be written as 



/dk ^ f d^k I P 

— — ^e-*'^^ (ofc cos ujkt + bk sin ujkt)+ / rfV— — e^'^^^"^ ) ^ J{x') 

The contribution of the principal value is: 

(^^0 ^iko(t-t') P ^ /" ^gifeo(i-t')_L ( ^ ^ 

27r -kl + ul J 27r 2uk\ko + uJk ko - c 

= ^S{t - t') sinujk{t - t') (C.l) 

where we have used 

dk- e'^" = / dk—^ — = e^^^i / dk—^^ = e'^^'inSix) 

k — a J k — a J x 

and S{x) = 6{x) — 9{—x) is the sign function. All this leads to 

dk 



[x] = I , ""y,^^ e *'''' (afc cos w^ t + 6^ sin cu^ t) + 

d-x' J j^^-^^ -^^^ S{t - t') sin cok {t - t') J{x') (C.2) 

In momentum space, 

0_fe(t) = ak cos cukt + bk sin Ukt + - — / d"'x' e^^^ S{t — t') sin ook {t — t') J{x') 

Swfc J 
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C.l Dirichlet boundary conditions 

The solution that vanishes at t = tj as well as at t = tf then reads 

dk dx' e~^''^^~^'^ 



{X 



dt'J{t', x) sin Uk (ti — t') S (t' — U) sin ojk (tf — t) + 

(C.3) 



(27r)"-i2a;fc sin LUk{tf - U 
+ sin ojk (tf — t') S (t' — tf) sin ojk (t — ti) + sin ujk(tf — ti) sin ojk (t — t') S (t — t' 



It vanishes for J = 0, in agreement with previous results. 

This means that the correct propagator to be used in the integral over V^ is 
given by 



V(x,x') 



dke "^^^^ ^'^ 



sin Uk {ti — t') S (t' — ti) sin Uk {tf — t) + 



2(271)" ^Uk sin WfcT 
sin ojk {tf — t') S {t' — tf) sin Uk {t — ti) + 

sin UkT sin Uk {t — t') S {t — t') 

That is, this is the only solution to the equation 

(n + m^) V{x, x') = 5{x- x') 

such that 

1^{x,x%^u = V{x,x%,^^^=0 

C.2 Feynman boundary conditions 

The boundary conditions for the Feynman propagator are defined by 

i(pk{tf) = ujk(t)k{tf) 
i<Pk{ti) = -u)k(pk{ti) 



(C.4) 



(C.5) 



In momentum space 






ikx 



4>-k{t) = — Wfcflfcsin CO kt + bkiOk COS cokt + / d"'x' S{t — t')i0kcos uJk{t — t')J{x') 

Jr" ^^k 

(the delta function does not contribute). 



The boundary conditions are then 



Akx' 



Uk ak cos Uk tf + bk sin Uk tf + 



2Wfc Jjjn 



rx'l 

rfV e*^"" S{tf - t') sinuk {tf - f) J{x') 



— iWfcOfcSin Uktf + ibkUkCos Uktf + i / d^x' - — S{tf — t')ukCos Uk{tf — t')J{x') 

Jm." ^^k 

1 
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Akx' 



—iukaksin Ukti + ibkCOk cos Ukti + i I dTx - — S(ti —t')ukCos Ukiti — t')J{x') 
—cuk I ttk COS Uk ti + bk sin Uk U H / d'^x' e^^^ S{ti — t') sin Uk {ti — t') J{x') 



This can be written as 

and this yields 

afc = ^ e-^"^^ /" d"a;'e^'^"" J(a;') (^(t/ - t') e— ^(*'-^) - 5(t, - t') e*-^(^+*') 
ci"a;'e^^^' J(a;') (^(t/ - t') e'^^^^' - S{U - t') e*^'=*' 



4wt 



(C.6) 



hk 



Aco. 



jn / ikx' T( rf'\„—ii^kT 



k JR" 



dTx'e"''' J(xOe-*"'=' S{tf - t') e*"'=^' "* ^ + S{ti - t') e 



/N iu)k{T~t') 



./^ ^Jc^fc(T+t') 



4Wfc 



dVe*^"" J(a:') {s{tf - f) e-''^"'' + S{ti - t') e'"'"'' 



(C.7) 



This means that the corresponding classical solution reads 

(Tx' 1^t{x,x')J{x') 



\x 



with Ay the finite time Feynman propagator: 



dk 



^t{x^-)^ lj^^^.--''''-'''^T{k) 



Arik) 



Auk 



S{tf - f) e-'^^'' - S{U - t') 6*^'=*' ) cos ookt+ 



+ i (S{tf - t') 6"^^'=*' + S{U - t') 6*'^'=*') sin Ukt 



+ 



+ - — S{t — t') sin ujk{t — t') 
2uk 



(Ci 



It is then plain that in the limit tf = —ti = '^ and T — )■ oo, Feynman's continuum 
result is recovered. 

lim Ar(A;) = Ai?(A;) = [cos uJk{t — t') — iS{t — t') sin ujk{t — t')] 
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In general it yields 



At(A;) 



4wi 



[S{tf-t') - S{ti - t') - 2] cosujk{t - t') + 

+ i[S(tf - t') + S{U - t')] fiinujkit - t')) + A^(A;) (C.9) 



The above results are valid for general sources with arbitrary support. When (as 
in our case) the support is restricted to the interval ti < t' < tf, it is quite easy to 
check that the finite time propagator with Feynman's boundary conditions coincides 
exactly with ine usual Feynman's propagator. 

C.3 Feynman's propagator including de Sitter boundary terms 

The boundary terns that appear when redefining the field in de Sitter space imply 
(after the splitting (p = (p^ + x) ^''^ addition of 

n- 2 n- 2 

^^^Mzf)x{zf) - i^^<Pc{zi)x{Z'i) = i-ffm(f)c{zf)x{zf) - i-i,^(p^{zi)x{zi) 

to the boundary. In order to eliminate those cross-terms the boundary conditions to 
be imposed are 



i^k{zf)-uJkl<Pl{zf)+i-ifml(pl{zj) = 
i4)l{zi) + Ukl <fk{zi) + i'^iml (plizi) = 



(CIO) 



Let us now can make an slightly different antsatz for the form of the propagator, 
namely 

Ax(A;) = S{z — z') sinukl{z — z') + 

2ujk 

+ Cfc cosUkl{z — z') + bk sinuJkl{z — z') (C.ll) 

The previous equations for the boundary values give rise to the following coefficients: 



-ilu:^Z 



D 



ook {-mim^ (e2i/a;feZ _ 1) _ 2imuJk{lf - li) + 4wD 
ak=^ {{2uk - 27/m)e2*'-'=(^/-^') + i^fm) [jm (-1 + e^^'^'^^^-^')) + 2tUk) 
■{S{zf - z') - S{z, - z'))e'^^k{2z'-zf-z^) 



Dk 

4 



T ^ -^k ^-il^^^{2z'+Zf+z,) 



+'^ime 



S{zf - z') (7/me=^''^'"'= + e'*'"'=^^(7/(-m) - 2iUk)) (e^''''^'{-i,m - 2iUk) + 

(C.12) 



;2*'^'"H7i^ - 2«Wfc) - 7»me2*'^'="" 
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where Z = Zj — Zi. 

However, the interaction takes place only in the interval \zi^Zj\^ so we should 
restrict the variables z and z' to be into this interval, so: 






Also, since the exponent of the Feynman Kernel is symmetric in the source J, 
we must symmetrise the propagator. 
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